A note on estimating stochastic volatility and its volatility: a new 

simple method 



With the exception of Alghalith (2012), the previous methods of esti- 
mating (stochastic) volatility suffered two major limitations. Firstly, they 
are only applicable to time-series data and therefore the volatility of the 
investor's portfolio (a cross section of assets) for time t cannot be properly 
and directly estimated. Secondly, these methods required the arbitrary gen- 
eration of data series for volatility in order to estimate the volatility. This 
seems somewhat contradictory and futile. Alghalith (2012) used a stochas- 
tic factor model, where volatility is a function of an external economic factor 
(GDP). In this note, we use a stochastic volatility (in the sense of Heston, 
etc.). In doing so, we overcome the aforementioned limitations. Below is a 
description of the theoretical model. 

We modify the standard portfolio model to include a stochastic volatility. 
We still use a two-dimensional Brownian motion {Wu, W2s, J^s}t<s<T 
probability space {Q,J^s, P) , where {J^s}t<s<T augmentation of filtra- 

tion. We include a risky asset and a risk- free asset. The risk-free asset price 

T 

process is given by 50 = 6* , where r E Cl (R) is the rate of return. 



The dynamics of the risky asset price are given by 



dSs = S, {njs + asdWis} , (1) 

where /i^ and as are the rate of return and the volatihty, respectively. 
As in Heston's model, the stochastic volatihty is given by 

dal = (a - /3(7^) ds + -fasdW2s: ol = a, (2) 

where |p| < 1 is the correlation factor between the two Brownian motions, 
a, 13 and 7 are constants. 

The wealth process is given by 

= a; + j {r,X; + (/x, - r,) ttJ ds + j T^sf^sdWu, (3) 
t t 

where x is the initial wealth, {t^s, J's}t<s<T portfolio process, and 

T 

E J TT^ds < 00. The trading strategy tt^ G A{x,a) is admissible. Xg = 
t 

TTg + Bs, where Bg is the amount invested in the risk- free asset. 

The investor's objective is to maximize the expected utility of the terminal 
wealth 



2 



V {t, X, a) =supE [U (XJ) I J^t] , 

TT 

where V {.) is the value function, U (.) is continuous, smooth, bounded and 
strictly concave utility function. 

The value function satisfies the Hamilton- Jacobi-Bellman PDE (suppress- 
ing the notations) 



Sup I ^7rVl4^ + [k{h- r)] ^ipai^V^^ \ = 0, 



V (T, x,a)^u {x, a) . 

Thus the solution yields 



TT = 377 PlTT- 

^ y TT * TT 



We take a second-order Taylor's expansion of V {t, x, a) and therefore 



Vx {t, X, cr) ~ ckq + OiiX + a2(T. 



Substituting this into (jlj yields 



^ - r) [ao + ai(7r* + S) + 020"] 02 
TT = z P7 — • 



The above equation can be rewritten as 

- r) TT* (/i - r) ao + ftiS + 02 

TT =P2 = \ = ;Pl = ,P2 = -P7 — , 

a a «! 

(5) 

where is a constant. Thus 



* , /3i(/i-r) 

^ = :; 71^ ^ • 

1 + ^ a + n-r 



The above equation can be rewritten as the following regression equations 



^2 +Mi^, (6) 



where /3j is a parameter to be estimated by a non-linear regression (while 
the variables vr*, /i and r are observed data), and is an estimate of the 
volatility of the portfolio for period t. 

To estimate the volatility of volatility 7, we multiply (E]) by 7 and take 
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the inverse to obtain 



where is an estimate of 7, $2 estimated volatihty from the previous 

regression, and tt is observed data. A similar procedure can be used to 
estimate the factor of correlation. 
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